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The effect of external electric field on electron-hole correlation in GaAs quantum dots is investigated. The 
electron-hole Schrodinger equation in the presence of external electric field is solved using explicitly corre- 
lated full configuration interaction (XCFCI) method and accurate exciton binding energy and electron-hole 
recombination probability are obtained. The effect of the electric field was included in the 1-particle single 
component basis functions by performing variational polaron transformation. The quality of the wavcfunc- 
tion at small inter-particle distances was improved by using Gaussian-type geminal function that depended 
explicitly on the electron-hole separation distance. The parameters of the explicitly correlated function were 
determined variationally at each field strength. The scaling of total exciton energy, exciton binding energy, 
and electron-hole recombination probability with respect to the strength of the electric field was investigated. 
It was found that a 500 kV/cm change in electric field reduces the binding energy and recombination prob- 
ability by a factor of 2.6 and 166, respectively. The results show that the eh-recombination probability is 
affected much more strongly by the electric field than the exciton binding energy. Analysis using the polaron- 
transformed basis indicate that the exciton binding should asymptotically vanish in the limit of large field 
strength. 
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I. INTRODUCTION 

The influence of external electric field on optical prop- 
erties of semiconductors has been studied extensively us- 
ing both experimental and theoretical techniques. In 
bulk semiconductors the shift in the optical absorbing 
due to the external field is known as the Franz-Kcldysh 
effect. 1 In quantum wells and quantum dots, applica- 
tion of electric field has shown to modify the optical 
properties of nanosystems and is known as the quantum- 
confined Stark effect (QCSE). 2, ! The application of the 
external field induces various modifications in the opti- 
cal properties of the nanomaterial including, absorption 
coefficient, spectral weight of transitions, and change in 
Amax of the absorption spectra. In certain cases, the ap- 
plied field can lead to exciton ionisation. 1 The quantum- 
confined Stark effect has found application in the field of 
electro-absorption modulators, solar cells 1 ' and the light- 
emitting devices. Recent experiments by Weiss et al. on 
semiconductor quantum dots have shown that the QCSE 
can also be enhanced by the presence of heterojunctions. 
In some cases, the QCSE can be induced chemically be- 
cause of close proximity to ligands. The QCSE also plays 
a major role in electric field dependent photoconduc- 
tivity in CdS nanowires and nanobelts/' Electric field 
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has emerged as one of the tools to control and cus- 
tomize quantum dots as novel light sources. In a re- 
cent study, electric field was used in generation and con- 
trol of polarization-entangled photons using GaAs quan- 
tum dots. It has been shown that the coupling be- 
tween stacked quantum dots can be modified using elec- 
tric field. 11 

The QCSE has been investigated using various the- 
oretical techniques including perturbation theory, 12-1 ' 
variational techniques, 18-24 and configuration interaction 
method. 20 In the present work, development of ex- 
plicitly correlated full configuration interaction (XCFCI) 
method is presented for investigating effect of external 
electric field on quantum dots and wells. The XCFCI 
method is a variational method in which the conven- 
tional CI wavefunction is augmented by explicitly corre- 
lated Gaussian-type geminal functions. !( ' The inclusion 
of explicitly correlated function in the form of the wave- 
function is important for the following two reasons. First, 
the addition of the geminal function increases the conver- 
gence of the FCI energy with respect to the size of the un- 
derlying 1-particle basis set.' 3 ' Second, inclusion of explic- 
itly correlated function improves the form of the electron- 
hole wavefunction at small inter-particle distances which 
is important for accurate calculation of electron-hole re- 
combination probability. 38-40 The effect of explicitly cor- 
related function on the convergence of CI energy has been 
investigated by Prendergast et al. and is directly re- 
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lated to accurate treatment of the Coulomb singularity 
in the Hamiltonian. 57,41 ' 42 Varganov et al. have demon- 
strated the applicability of geminal augmented multi- 
configuration self-consistent field wavefunction for many- 
electron systems. 1 ' Elward et al. have also performed 
variational calculation using explicitly correlated wave- 
function for treating electron-hole correlated in quantum 
dots. 40,44 

One of the important features of the XCFCI method 
presented here is the inclusion of the external field in the 
ansatz of the wavefunction. This is achieved by defining 
a new set of field-dependent coordinates which are gener- 
ated by performing variational polaron transformation 
and recasting the original Hamiltonian in terms of the 
field-dependent coordinates. The variational polaron 
transformation was introduced by Harris and Silbey for 
studying quantum dissipation phenomenon in the spin- 
boson system 4 '' and is used in the present work because 
of the mathematical similarity between the spin-boson 
and the field-dependent electron-hole Hamiltonian. 

The remainder of this article is organized as follows. 
The important features of the XCFCI method are sum- 
marized in Sec. II A, construction of the field dependent 
basis functions is presented in Sec. II B, the application 
of the XCFCI method using field-dependent basis is pre- 
sented in Sec. Ill, and the conclusion are provided in 
Sec. IV. 



II. THEORY 

A. Explicitly correlated full configuration interaction 

The field dependent electron-hole Hamiltonian is de- 
fined as 46 ' 4 ' 



H = 



2m Q 
1 

e|r e h| 



c 2m h 11 
|e|F-(r e -r h ) 



(1) 



where m G is the mass of the electron, is the mass 
of the hole, e is the dielectric constant, and F is the 
external electric field. The external potential v° xt and 
u^ xt represent the confining potential experienced by the 
quasi-particles. The form of the XCFCI wavefunction is 
defined as 



*xcfci = G^c fc $ A 



(2) 



where is the CI coefficient and are basis functions. 
The operator G is known as the geminal operator and is 
an explicit function of r c h and is defined as 



N c N h N g 
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where N g is the number of Gaussian functions included 
in the expansion, N e and N c are the number of electrons 
and holes, respectively. The parameters and 7^ used in 
the definition of the geminal operator are obtained vari- 
ationally. The construction of the basis functions used 
in the definition of XCFCI wavefunction in Eq. (2) will 
be discussed in Sec. II B. The XCFCI calculation is per- 
formed in two steps. In the first step, the parameters of 
geminal operator are obtained variationally by perform- 
ing the following minimization 



. (G$ |g|G$ ) 
E[G ^ ] = (G*o|G*o) 



(4) 



In the second step, the expansion coefficients {c^} are 
obtained variationally and are defined by the following 
minimization procedure 
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The above equation can be rewritten as a FCI calculation 
of transformed operators 
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where the transformed operators are defined as 
1 
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The exact expression of the transformed operators in Eq. 
(7) and (8) and discussion relevant to their derivation 
has been presented earlier in Ref. 40 and 44 and is not 
repeated here. The -Excfci reduces to conventional FCI 
energy in the limit of geminal function equals to 1 



Etci — Jim -Excfci 



(9) 



We expect the -Excfci energy to be lower than the FCI 
energy for identical set of basis functions and earlier stud- 
ies have shown this to be true. 44 

After the successful completion of the XCFCI calcu- 
lations, the field dependent exciton binding was calcu- 
lated from the difference between the non-interacting 
and interacting ground state energies. Defining the non- 
interacting Hamiltonian as 



H = lim H, 

e -1 -)-0 

the exciton binding energy is computed as 

,(0) 
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12) 

¥|flo|*) 



where Eq is defined in Eq. (12) 
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The field dependent electron-hole recombination proba- 
bility is obtained from the XCFCI wavefunction using the 
following expression 40 ' 44 

Pi [F] = ^gggil^ ~ r h)I^XCFCl) , 13 > 
(^XCFCll^XCFCl) 

The exciton binding energy and the recombination prob- 
ability arc functionals of the applied external field and are 
indicated explicitly in Eq. (11) and (13), respectively. 

B. Construction of field dependent basis set 

One of the key features of the electron-hole Hamilto- 
nian used in the present work is the presence of the field- 
dependent term in Eq. (1). Since the convergence of the 
CI expansion depends on the quality of the underlying 1- 
particle basis, it is desirable to construct and use efficient 
single particle basis sets. In the present work, we have de- 
veloped field-dependent basis functions and the details of 
the derivation are presented as following. Starting with 
the expression of Ho in Eq. (10), the zeroth-order Hamil- 
tonian is expressed as a sum of non-interacting electronic 
and hole Hamiltonians 

H = H c a + ff h , (14) 

where the expression for the single-component non- 
interacting Hamiltonian is given as 

H c o=T c + v^ + |e|F-r c (15) 
H% = T h + v^-\e\F-r h . (16) 

As seen from the above equation, the coupling between 
the external field and the quasiparticle coordinates is lin- 
ear. The above Hamiltonian shares mathematical sim- 
ilarity with the spin-boson Hamiltonian that has been 
used extensively in quantum dissipative systems. In the 
present method, we perform analogous transformation 
which is defined by the follow equations 

q c = r c + A C F (17) 
q h = r h - A h F. (18) 

Similar to the polaron transformation in the spin-boson 
system, the coordinates of the quasiparticle experience 
a shift due to the presence of the external field. 48 Us- 
ing the method of variational polaron transformation by 
Harris and Silbey, the shift parameter A is determined 
variationally. The field-dependent electronic basis func- 
tions are obtained by first constructing the Hamiltonian 
matrix using Gaussian-type orbitals (GTO) and then di- 
agonalizing the resulting matrix 

H^ = e?(A e )$f i=l,...,M e (19) 
H^ = ^(X h )^ j 1 W„. (20) 



The value of the shift parameter is obtained variationally 
by minimizing the trace 

min^ef A c . (21) 

i 

The Ah is also obtained by a similar procedure. The 
electron-hole basis functions for the FCI calculations are 
constructed by taking a direct product between the set 
of electronic and hole single-component basis sets 

{$ fc } = {$°} g, (22) 

The procedure described above is a general method that 
is independent of the exact form of the external poten- 
tial. However if the external potential is of quadratic 
form, the field dependent zeroth-order single-component 
Hamiltonian has an uncomplicated mathematical form 
and additional simplification can be achieved. 

III. RESULTS AND DISCUSSION 

The electron- hole Hamiltonian in Eq. (1) 
has been used extensively for studying optical 
rectification 15 ' 1 '' 46 ' 47 ' 49 effect in GaAs quantum dots 
and all the system specific parameters were obtained 
from previous calculations on the GaAs system. 46 ' 47 The 
parabolic confinement potential has found widespread 
applications 12,15 ' 20 ' 2 ' 3 ' 50-61 in the study of quantum dots 
and was used in the present work to approximate the 
external potential term in the Hamiltonian. All the 
parameters that are needed for complete description of 
the electron-hole Hamiltonian used in the calculations 
are presented in Table I. Following earlier work on the 

TABLE I. System dependent parameters used in the electron- 
hole Hamiltonian for the GaAs quantum dot ' 



Parameter 


Value 


m 


0.067m 


TOh 


0.090m 




9.048 x 10" 7 a.u. 


fch 


1.122 x 1CT 6 a.u. 


C 


13.1e 



effect of electric field on non-linear optical properties 
of GaAs quantum dots, ' the external electric field 
was aligned along the z-axis and the field strength 
was varied from zero to 500 kV/cm. Similar to the 
spin-boson Hamiltonian, the polaron transform resulted 
in shifted harmonic oscillators. The eigenvalues and 
and eigenfunctions of the Ho were obtained analytically, 
and the lowest ten eigenstates of the shifted harmonic 
oscillator Hamiltonian were used in the construction 
of the 1-particle basis. The direct product between 
the electronic and the hole basis sets was performed to 
generate the electron-hole basis for the FCI calculations. 
The geminal minimization was performed using a set of 
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three {6fc,7fc} parameters at each field strength, and the 
optimized values are presented in Table II. The total 
exciton energy for the field-free case was found to be 
269.45 meV. The total exciton energy of the system as 
a function of the field strength is presented in Fig. 1 . It 
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FIG. 1. Relative exciton energy compared to the fit E = 
(-2.7925 x 10~ 6 )F Z 2 + (-7.0938 x 10~ 5 )F Z + 1. 

is seen that the total energy decreases with increasing 
field strength. Earlier studies on this system indicate 
that the exciton energy is a quadratic function of the 
applied field. 62 ' 63 To investigate the scaling of the total 
exciton energy with respect to the field strength, we have 
performed least-square fit of the calculated values with a 
second order polynomial and the results are presented in 
Fig. 1. The results from these calculations confirm that 
the quadratic scaling of the exciton energy as a function 
of the field strength. The exciton binding energy was 
calculated using Eq. (11) and was found to be 28.52 
meV for the field-free case. The effect of the external 
field on the exciton binding energy was investigated by 
calculating the relative binding energy which is defined 
by the following equation 



E B 



E b [F] 
E B [F = 0Y 



(23) 



It is seen from Fig. 2 that the exciton binding energy 
decreases with increasing field strength. As the field 
strength is increased from to 500 kV/cm, the exciton 
binding energy decreases by a factor of 2.6. In addition 
to calculation of binding energy, the effect of the field 
on electron-hole recombination probability was also in- 
vestigated. Analogous to the relative binding energy, the 
relative recombination probability is defined as 



P. 



eh 



Peh[P = 0]' 



(24) 



and is presented in the Fig. 2. It is seen that the there is a 
sharp decrease in the recombination probability with in- 
creasing field strength and the recombination probability 
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FIG. 2. Comparison of Eb and P e h as a function of electric 
field strength. 



at 500 kV/cm is lower than the field-free case by a factor 
of 166. One of the key results from this study is that exci- 
ton binding energy and eh-recombination probability are 
affected differently by the external electric field. It is seen 
that the exciton binding energy and eh-recombination 
probability follow different scaling with respect to field 
strength. 

The polaron transformation also provides insight into 
the effect electric field on the exciton binding energy 
in the limit of high field strengths. Starting with the 
transformation defined in Eq. (17), the electron- hole 
Coulomb interaction in the transformed coordinate can 
be expressed as 



1 



(q c - «&) - (A e + X h )F\ 



woh(q). (25) 



It is seen in the above equation that the above expres- 
sion will be dominated by the field-dependent term in the 
limit of high field strength. A direct consequence of the 
above condition is that in the limit of high field strengths, 
we expect the exciton binding energy to be small 



ff(q) « ffo(q) 



Et- 



fori < |F| < oo. (26) 



It is important to note that the above conclusion is inde- 
pendent of the choice of the external potential. 



IV. CONCLUSION 

The effect of external electric field on exciton bind- 
ing energy and electron-hole recombination probability 
was computed using explicitly correlated full configura- 
tion interaction method. Field-dependent basis functions 
were used in the calculations and a variational polaron 
transformation scheme was developed for construction of 
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TABLE II. Optimized geminal parameters used in the calculations of energy and recombination probability. 

F z (kV/cm) 100 200 300 400 500~ 



bi 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


71 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


b 2 


1.40 x 10 _1 


9.99 x 10" 1 


4.99 x 10~ 2 


5.78 x 10~ 3 


1.00 x 10~ 2 


5.59 x 10~ 3 


72 


2.29 x 10~ 4 


4.60 x 10~ 6 


1.11 x 10~ 2 


1.11 


1.00 


1.11 


63 


4.35 x 10~ 2 


1.08 x 10 _1 


8.90 x 10~ 2 


1.67 x 10~ 2 


2.00 x 10~ 2 


1.58 x 10~ 2 


73 


1.13 x 10~ 2 


1.00 x 10~ 2 


1.01 x 10~ 3 


1.11 x 10 _1 


1.01 x 10 _1 


1.02 x 10 _1 



field-dependent basis functions. It was found that both 
exciton binding energy and electron-hole recombination 
probability decrease with increasing field strength. One 
interesting conclusion from this study is that the bind- 
ing energy and recombination probability follow different 
scaling with respect to the external electric field. For the 
range of field strengths studied, the recombination prob- 
ability and exciton binding energy decrease by a factor of 
166 and 2.6, respectively. These results give important 
insight into the application of electric field for manipu- 
lating excitons in quantum dots. 
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